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Abstract 

Tensor calculation of suffix-contraction is carried out by a C-program. Ten- 
sors are represented graphically, and the algorithm makes use of the topology 
of the graphs. Classical and quantum gravity, in the weak-field perturbative 
approach, is a special interest. Examples of the leading order calculation of 
some general invariants such as R^ u \ (7 R fJ ' u> ' cr are given. Application to Weyl 
anomaly calculation is commented. 
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1 Introduction 



The field theory is described by a Lagrangian which has fields of different spins and 
their different derivatives: ip, d^ip, d^duip, A^, d u A^, h^, d a h^ V) ■ ■ ■ (<p : scalar, 
A^ : vector, h^ u : (2nd rank) tensor ) . For the n- dimensional flat Euclidean(Lorentzian) 
space, a general form is a k-th rank tensor ( we call it k-tensor in the following, 
k=0,l,- • •) under the global SO(n) (SO(n-l,l)) symmetry. In the electro-magnetism, 
the energy density is given by, in the Euclidean space, an SO(n)-invariant: 

E 2 + B 2 = -F^F^v 

= (d,A„) 2 - d,A u ■ d„A, , (1) 
= dpAjj — d u A^ 

Similar thing can be said for the Yang-Mills theory. For the curved space, it is 
valid when we treat a gravitational theory in the weak-field perturbation : g^ v = 
8fj,v + h^ v , \hp V \ <C 1. The scalar curvature R is given as 

R = d ll d fl h uu -d fl d l/ h flu + 0(h 2 ) . (2) 

This is a 0-tensor. Generally physical quantities are expressed by invariants (0- 
tensors) and all suffixes (both derivative- suffixes and field-suffixes) are contracted 
. It is sometimes, however, cumbersome to do calculation involving the suffix- 
contraction, especially for a higher-spin fields such as gravitational ones. 

The importance of high-computational ability in the development of the quan- 
tum gravity is historically shown as follows. In 1963 Feynman|l[] calculated some 
scattering amplitudes involving gravitons (, which turned out to be an important 
finding in gauge theories: Feynman-DeWitt- Faddeev-Popov ghost) by doing alge- 
braic computer calculation. In 1974, 'tHooft and Veltmanp| fully computed the 
1-loop ( both matter-field loop and graviton loop) counter-terms of gravitational 
theories for the first time (by 'finger' calculation using a highly-refined calculational 
method called 'background- field method'). This result means the quantum gravity 
is (1-loop) unrenormalizable when matter fields are coupled. It urged theoretical 
physicists, in the latter half of 70's, to look for or examine other alternative grav- 
itational models such as supergravity. In 1985, Goroff and Sagnotti|| succeeded 
in calculating the full 2-loop counter-terms of the pure gravity using a C-language 
program. Their result was re-confirmed by van de VenQ in 1992 by use of a nice 
algebraic software called 'FORM'||. The gauge dependence of their results was ex- 
amined by H using a REDUCE program. The result that 2-loop counter-terms do 
not vanish on-shell let many people have an opinion that Einstein gravity could not 
be regarded as the fundamental quantum theory. The present high activity in the 
string theory is based on these past experiences. The above history of these three 
and half decays clearly shows, for important development of gravity, good algebraic 
computer programs and efficient computational techniques play important roles. 

2 Generally the summation over all space-suffixes which makes an inner product, such as A^A^ = 
52u=i A^A^, is called Einstein's contraction. 
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This paper is partly motivated by the analysis of the Weyl anomaly. Compared 
with the chiral anomaly, the general structure of the Weyl anomaly is not so simple. ( 
For the chiral one, there is a result for the general n dimension [0.) There exist some 
analysis || H |K], |T| and the most popular belief is that they are given by, at each 



(even) dimension, the Euler term, some Weyl invariants and some 'trivial terms'. 
The explicit results are obtained only up to the 4 dimensional case ||12fl . Clearly 
a way to improve the situation is to explicitly obtain the Weyl anomaly in higher 
dimensions. We tried it at the 6 dim, but we soon realize its calculation involves so 
much tensor contraction that is difficult in the hand- calculation. Furthermore, the 
popular algebraic softwares are inconvenient or inefficient for the problem (] . 

In this paper, we present a new algorithm for the tensor contraction and present a 
C-program ( we call it "WEAKGRAV") using it. Important tensors often have high 
permutation symmetries with respect to their suffixes. For example, the Riemann 
tensor has 

R^iuafi RaP^iu RufiafB Rfiufia • (3) 

When we treat 'products' of these tensors, it is most efficient to express the tensors 
by graphs which have their permutation symmetries |i~3"H . All tensor products are 



completely specified by the topology of graphs. It makes one free from bothering 
about suffixes. We have only to deal with topology of graphs, not their suffixes. The 
present algorithm makes use of this idea. In HT3| , the general tensors in the curved 
space are treated. In the present paper, we treat global SO(n) tensors in the flat 
space. 

Generally the local properties of the quantum gravity can be analysed by the 
perturbation around the flat space. For example, the Weyl anomaly , the effective 
action and the counterterms are obtained. ( We use, in the final form, the fact 
that they must be invariant under the general-coordinate transformation. This 
procedure allows us to take into account all orders of the weak-field expansion.) 
In the case, there appear many terms of "products" of global SO(n) tensors. The 
present algorithm is useful in such cases. 

The main purpose of this paper is to explain how the above idea about tensor- 
contraction is realized in the explicit computer calculation. We give two sample 
calculations. Emphasis is not on the calculated results themselves but on how the 
results are obtained. Such presentation is good for readers to apply the present 
algorithm to their own problems. 

Finally we describe the main features of the program "WEAKGRAV" . 

• Computers: all computers running C-program 

• Program language used: C-language 

• No. of lines including sample data: 300 

• Input files needed: indata.dat 



3 Most softwares so far treat suffixes directly and the permutation symmetry w.r.t. suffixes 
is coded by explicitly writing (or declaring) the equality between quantities with the permuted 
suffixes. This algorithm is not appropriate in the higher dimensions or in higher orders. See Sec. 6. 
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Figure 1: 4-tensor d^d^hap 

• Purpose of program: In the perturbative treatment of gravitational theories, we 
face many terms with complicated suffix contraction. This program helps to do 
such calculation very fast. Each term is graphically represented. Topological 
indices of graphs are utilized in this algorithm. It can be applied to general 
field theories. 

• Restrictions on the complexity of the problem: The program is at a basic level. 
The contraction terms are restricted to the quadratic type ones: dfj,d u h\ a ■ 
d a dphys- This is, however, sufficient to do Weyl anomaly calculation of 4 dim 
quantum gravity. 

The program "WEAKGRAV" is available on the request to the author's home in- 
stitute. ( It is restricted only to non-commercial use.) 

The content is organized in the following way. In Sec. 2 we explain the graphical 
representation of global SO(n) tensors, which is used in the explanation of the 
algorithm. The codes for the tensors are explained in Sec. 3. In Sec. 4, new quantities 
called ''indices' 1 are introduced and explained. They express the topology of global 
SO(n) invariants. We demonstrate the calculation using two simple examples in 
Sec. 5. We conclude in Sec. 6. In App.A, all (<9<9/i) 2 -invariants are graphically shown. 
Finally, in App.B, we add the input and output data for the sample calculations of 
Sec. 5. 

2 Graphical Representation of Global SO(n) Ten- 
sors 

In the programming explanation, we use a graphical representation for a tensor 
because it can clearly express the connectivity of suffixes in a tensor |T^, |TJJ]. In 
this section we explain the minimum amount of the representation necessary for the 
present programming. 

We can graphically express a basic element of the present program: 4-tensor 
d^d u h a p as in Fig.l. We call the dotted lines suffix-lines, the rigid line bond and 
the intersections of two dotted lines and one bond vertices. There are two kinds of 
vertices: the vertex with a crossing symbol 'x' and that one without it. The suffixes 
of derivatives (/i, v in Fig.l) are associated with those suffix-lines which extend from 
a vertex without the crossing symbol, whereas the suffixes of fields (ex, [3 in Fig.l) 
are associated with those which extend from a vertex with it. This graph respects 




Figure 2: 2-tensors of (a) d 2 h a p , (6) d^d u h aa and (c) d^dph^ 




Figure 3: Invariants of P = d^d^haa and Q = d a dph a p . 

all symmetries of d^d u h a p w.r.t. the permutation of suffixes: d^d v h a p = dyd^hap = 
d^duhpa. The suffix contraction is expressed by connecting the same suffix-lines. 
For example, 2-tensors : d 2 h a p , d^d u h aa , d^dph a p , which are made from Fig.l 
by connecting two suffix-lines, are expressed as in Fig. 2. 

Two independent invariants (0-tensors) : P = d^d^h aa , Q = d a dph a p , which 
are made from Fig. 2 by connecting the remaining two suffix-lines, are expressed as 
in Fig.3. They are all possible invariants of ddh-type. Generally all suffix-lines of 
invariants are closed. Those closed suffix-lines are called suffix-loops. 

As for those tensors or invariants which are made from two or more <9<9/i-tensors, 
their graphical expressions are obtained in the same way. (See the following sec- 
tions.) We see a tensor and the corresponding graph are one-to-one. 



3 Code for Graphs 
3.1 Code for ddh- Tensors 

The basic element is the 4-tensor : d^h^. Here and in the following we change the 
Greek symbols, /i, v, a, /3, etc., to the non-negative integer numbers, 0, 1, 2, 3, 4, etc., 
for the convenience of programming. The four suffixes, characterizing this tensor, 
are stored in a 3 dimensional array of size 10 x 2 x 2 : d2h[i][j][k] as 

d2h[i][0][0] = 1 
d2h[i][0][l] = 2 
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Figure 4: 4-tensor of d^hu. 
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Figure 5: Invariants of P = did\h 2 2 and Q = d^h^ . 

d2h[i][l][0] = 3 
d2h[i][l][l]=4 

where we regard this tensor as "i"-th 'bond' (<9<9/i-tensor), "j" specifies a chosen 
vertex-type (j=0 for the vertex without x, j=l for the one with x) of the i-th 
ddh-tensor and "k" specifies one of two suffix-lines at each vertex (i,j). | See 
Fig.4. We call "i" 5 "j" and "k" as bond number, vertex-type number and suffix number 
respectively. 

In the examples of Fig.5, we store two invariants, P = did\h 2 2 and Q = <9i<92^i2 
in i=0 and i=l component of d2h[i][ ] ] respectively. 

d2h[0][0][0] = 1 d2h[l][0][0] = 1 

d2h[0][0][l] = 1 d2h[l][0][l] = 2 

There exists arbitrariness in the choice of the input data due to the suffix-permutation sym- 
metry of the tensor d^duhap. Because the graph identification is programmed in terms of its 
topology, we may take any one, say, 

d2h[i][0][0] = 2 
d2h[i][0][l] = 1 
d2h[i][l][0] = 4 
d2h[i][l][l] =3. 

5 As the size of "i", we here take 10. This allows us to treat up to 10 different 99/i-tensors at 
once. 
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Figure 6: 8-tensor of dxd^h^ ■ d^d^h^. 

d2h[0][l][0] = 2 d2h[l][l][0] = 1 

d2h[0][l][l] = 2 d2h[l][l][l] = 2 

Generally, contracted suffixes are dummy ones, that is, any suffix-numbers may be 
taken unless the same number is used for different contractions. In both examples 
above, the suffixes 1 and 2 are dummy ones and they can be substituted by any 
non-negative integers, say, 3 and 4 respectively. 

3.2 Code for (<9<9/i) 2 -Tensors 

The most general (<9<9/i) 2 -type tensor is the 8-tensor: <9i<9 2 /i 34 - d 5 d 6 h 78 which is graph- 
ically shown in Fig. 6. The code for Fig. 6 is 

d2h[0][0][0] = 1 d2h[l][0][0] = 5 

d2h[0][0][l] = 2 d2h[l][0][l] = 6 

d2h[0][l][0] = 3 d2h[l][l][0] =7 

d2h[0][l][l] = 4 d2h[l][l][l] = 8 

A 4-tensor d^h^ ■ d^d^h^, which is made from the above 8-tensor (Fig. 6) by 
contracting(connecting) the suffixes (suffix-lines) 3 and 4 with 7 and 8 respectively, 
is represented as in Fig. 7. The code for Fig. 7 is 

d2h[0][0][0] = 1 d2h[l][0][0] = 5 

d2h[0][0][l] = 2 d2h[l][0][l] = 6 

d2h[0][l][0] = 3 d2h[l][l][0] = 3 

d2h[0][l][l] = 4 d2h[l][l][l] =4 

where two indices 3 and 4 are dummy. There are some other 4-tensors which can 
be obtained from Fig. 6 by connecting two pairs of suffix- lines in different ways. 

An invariant: d\dih^ • <9i<92^34, which is made from Fig. 7 by contracting suffixes 
1 and 2 with 5 and 6 respectively, is graphically obtained as in Fig.8. The code for 

6 Because the numbering of two bonds is arbitrary ( permutation symmetry of two 99/i-tensors 
in a (<9<9ft,) 2 -tensor ), we may take, say, 

d2h[0][0][0] = 5 d2h[l][0][0] = 1 

d2h[0][0][l] = 6 d2h[l][0][l] = 2 

d2h[0][l][0] = 7 d2h[l][l][0] = 3 

d2h[0][l][l] = 8 d2h[l][l][l] = 4 
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Fig. 8 is 
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Figure 8: Invariant of B3 = d^h^ ■ didzh^. 



d2h[0][0][0] = 1 d2h[l][0][0] = 2 

d2h[0][0][l] = 2 d2h[l][0][l] = 1 

d2h[0][l][0] = 4 d2h[l][l][0] = 3 

d2h[0][l][l] = 3 d2h[l][l][l] =4 

All possible and independent invariants of (<9<9/i) 2 -type are listed in Appendix A. 
As mentioned before, all the suffix-lines of the invariants are closed. 



4 Indices for Characterizing (<9<9/i) 2 -Invariants 

By virtue of the graphical representation, we can identify a tensor by the topology of 
the corresponding graph. We need not to know the explicit suffix-names associated 
with the tensor any more. In the computer calculation, this approach reduces the 
memory space and speeds up the tensor calculation considerably. Let us explain 
how to do it concretely using the codes explained in the previous section. 
The following two items are the important steps for the tensor calculation. 

• Listing up all independent invariants (closed graphs). 

• Characterizing every closed graph by a set of some appropriate indices, which 
are computed from the code of the graph. Those indices must be uniquely 
defined irrespective of arbitrariness of the code due to 1) the suffix permuta- 
tion symmetry, 2) arbitrariness in numbering bonds, suffix-loops, tadpoles(see 
later), etc, 3) arbitrariness of dummy suffixes, 4) others (see later). 
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Figure 9: loopstreamf [ ][ ] for 53. Arrows indicate the directions of tracing along 
each loop. 



The former is generally done by the graph theory ||15|| , while the latter can be done 



by a program using the present algorithm. The indices are explained in this section. 

Let us begin with the simplest case of <9<9/i-invariants. See Fig.5. We denote, 
[the number of suffix-loops] — 1, as loop in the program. f\ The index of loop is 
sufficient for discriminating two independents, P and Q; loop=l for P = dydxh-xz and 
loop=0 for Q = d\d2hi 2 . In this paper we focus only on (<9<9/i) ^invariants in order 
to explain the algorithm most simply. The generalization to the general (ddh) m - 
invariants (m=l,2,3,- • •) is straightforward. (See Sec. 6.) For the (<9d/i) 2 -invariants, 
we can easily check, without the knowledge of the graph theory, that 13 independent 
invariants listed in App.A exhaust all possible independent ones|b|, |16|. Because 
all suffix-lines for invariants are always closed, we seek some good indices associated 
with each closed loop. Of course the number of closed loops, loop+1, is one good 
index. It is, however, not sufficient to discriminate 13 invariants completely. 

4.1 loop, loopstream[ ][ ][ ] and loopno[ ][ ] 

A 3-dim array loopstream[maxloopno][2][2], defined below, is initialized as 

loopstream[/][i][j] = 99, 
1=0,1,- ■ ■ <maxloopno=4; i=0,l; j=0,l 

For Z-th loop, we assign an ascending natural number, 1,2,3,- ■ ■, at each vertex (i,j) 
along the loop: loopstream[Z] [i] [j] = l,2,3,- ■ -. When a vertex (i,j) belongs to Z-th 
loop, we assign as loopno[i][j]=Z. See Fig. 9 for example. 

loop=l 

loopstream[0][0][0] = 1 loopno[0] [0]=0 

loopstream[0][0][l] =99 loopno[0] [1] = 1 

7 In the following, when the variable names, used in the program, appear in ordinary sentences 
of the text, they are underlined (' '). 

8 The choice of the starting vertex is arbitrary. This is one of "4) others" in the list of arbi- 
trariness, in the beginning of this section. The number "i" refers to "i-th bond" and it runs from 
to 1 because we are considering a quadratic term of ddh. The number "j" refers to a vertex-type: 
j=0 means the vertex without the crossing symbol, j=l means the vertex with it. See Fig.l. 
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Bl 



Figure 10: Invariant of Bl. # of vertices of 0-th and 1-th loop are 1 and 3 
respectively. 0-th loop is a tadpole. 

loopstream [0] [ 1] [0] = 2 loopno [1] [0] =0 

loopstream[0][l][l] =99 loopnojl] [1] = 1 
loopstream[l][0][0] = 99 
loopstream[l][0][l] = 1 
loopstream[l][l][0] = 99 
loopstream[l][l][l] = 2 

4.2 maxver [ ], vertex [ ][ ][ ] 

maxver [I] stores the number of vertices which Z-th loop has. When maxver [Z'] = l, 
the I'-ih loop is called tadpole. See Fig. 10. maxver [ ] for Fig. 10 is given by 

maxver [0] = 1 
maxver [1] = 3 

For later convenience, the same information as loopstream[ ] | ] is stored in another 
form: vertex [ ] ] | ] . When we trace the 1-th loop in a certain direction [] , let the 
k-th vertex (k=0,l,- • ■) be (i,j). Then we assign as vertex [Z] [k] [0]=i, vertex [Z][k][l]=j, 
k=0,l,- • • < maxver [Z]. For Fig.9 we have 

vertex[0][0][0] = 

vertex[0][0][l] = 

vertex[0][l][0] = 1 

vertex[0][l][l] = 

vertex[l][0][0] = 

vertex[l][0][l] = 1 

vertex[l][l][0] = 1 

vertex[l][l][l] = 1 

4.3 absdeli [ ], absdelj[ ], tadpoleno and tad[ ][ ] 

The number of tadpole-loops in a graph is stored in tadpoleno. When the t-th 
(t=0,l,- • • <tadpoleno) tadpole-loop connects with a bond at a vertex (i,j), we 
assign tad[t][0]=i, tad [t][l]=j. For CI, C2 and C3 of Fig. 11, we have the following. 

9 The choice of the direction of tracing is another arbitrariness in "4) others" in the beginning 
of this section. 
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Figure 11: tadpoleno and tad [ ] ] for CI, C2 and C3. 

CI C2 C3 

tadpoleno 2 2 2 

tad[0][0] 

tad[0][l] 1 1 

tad[l][0] 1 1 1 

tad[l][l] 1 

Now we come to key quantities in this algorithm. When we trace a suffix-line, 
along a loop, starting from a vertex (io Jo) m a certain direction, we pass some 
vertices, (ii,ji),(i2j2)r • • an d finally come back to the starting vertex. We focus 
on the change of the bond number (i) and the vertex-type number (j) when we pass 
from a vertex to the next vertex in the above tracing. See Fig. 12. For Z-th loop, we 
assign as E a l ong /_l 00 p l Ai l = absdeli [Z], E a l ong Z-loop l A jl = absdelj[/]. (They are 
called bond changing number said vertex changing number respectively, in Ref.||17||.) 
As examples, the indices for Al, A2 and A3 of Fig. 13 are given as 

Al A2 A3 

loop 

absdeh[0] 4 2 2 

absdelj jp] 2 2 4 

absdeli [ ] and absdelj [ ] defined above satisfy the following important properties 
for programming. 

1. They do not depend on the starting vertex for tracing along a loop. 

2. They do not depend on the direction of the tracing. 
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Figure 12: Change of i (bond number) and j (vertex-type number). 
Arrows indicate directions of tracings. 
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Figure 13: Calculation of absdelif ] and absdeljf ] for A1,A2 and A3. (\Ai\, \Aj\) 
are listed. 
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4.4 Table of Indices for all invariants 

We give the full list of indices for all invariants in Table 1. The program, WEAK- 
GRAV, computes all indices for each term, and identify it with a graph following 
the table. 
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Graph, (Run. No) \ Indices 


loop 


tadpoleno 


tad[][l] 


absdeli[ ] 


absdeljf 


Al = dffdxh,,,, ■ d a d u h^M 








nothing 


4 


2 


A2 = d a d x h Xfl ■ dadyh^l) 








nothing 


2 


2 


A3 = d a d x h Xfl ■ d^d u h ua ,(2) 








nothing 


2 


4 


Bl = d v d x K a ■ 0x8^,(3) 


1 


1 


1 


/ 


/ 


B2 = d 2 h Xu -d x d,h^(4) 


1 


1 





/ 


/ 


B3 = d^d v hxa ■ d^hxaX 5 ) 


1 





nothing 


9 


n 

u 


o 


n 






BA = d^dvhxc ■ dxdah^fi) 


1 





nothing 


9 


9 


2 


2 






Q 2 = (d^h^) 2 ,(7) 


1 





nothing 


o 


2 





2 






CI = d^d v h xx ■ dnd v h atT ,(8) 


2 


2 


i 


n 

u 


n 

u 




9 


n 


i 


n 


n 


C2 = d 2 V • d 2 V,(9) 


2 


2 


n 


n 
u 


n 
u 




9 


n 


n 


n 
u 


n 


C3 = d,d u h xx -d 2 h fiu ,(10) 


2 


2 


1 

-L 


n 

w 


n 













2 


2 


PQ = d 2 h xx -d,d v h^,(ll) 


2 


2 


1 






















2 


P 2 = (d 2 hxx) 2 ,(12) 


3 


/ 


/ 


/ 


/ 



Table 1 List of indices for all (<9<9/i) 2 -invariants. The symbol '/' means 
'need not be calculated'. Graph names, Al, A2, ■ ■ ■ are defined in Appendix A. 
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Note here that the quantities in Table 1 do not depend on the arbitrariness cited 
in the beginning of this section. This is one of key points of the present algorithm. 

5 How to use WEAKGRAV 

Let us explain how to use WEAKGRAV taking some simple examples. The program 
reads all input data from a file (named 'indata.dat'). The input and output data 
for two examples below are provided in Appendix B. 

Ex.1 R\ ilva B X[1,vu 

The weak-field expansion (g M „ = 5^ v + h^ v , |/t M „| <C 1) of Riemann tensor is given 

2R Xllvo = 2R x ^ + 0'(h 2 ) 

= d^hxa - d\d u h^ - v <-> a + 0{h 2 ) . (4) 

Changing the Greek suffixes (fj,,u,X,a) to non-negative integers, say, (1,2,3,4), we 
can get the input data of (termnoa, weightd2ha[ ],d2ha[ J ][ ]) for R\^ ua and 
(termnob, weightd2hb[ ],d2hb[ ][][]) for R x » va . 
The output says 

2R Xilvu x 2R x ^ a = -8 x Al + 4 x B3 + 4 x BA + 0(h 2 ) , (5) 

where Al = d a d\h^ u -d a d u h^ B3 = d^d u h\ a -d^d u h\ a , BA = d^d u h\ a -d\d a h^ (Table 
1, Appendix A). 



Ex.2 R^RT 

R pv and R^ v is expanded as 

2R IXV = 2R^ + 0'(h 2 ) 

= d^d u h aa - d^d a h au - d v d a h afl + d 2 h liv + 0{h 2 ) . (6) 

Changing the Greek suffixes (/x, v, a) to (1,2,3) for termnoa, weightd2ha[ ] and d2ha[ 
|||] (input data for 2R^ U ) and changing (/i, u, a) to (1,2,4) for termnob, weightd2hb[ 
] and d2hb[ ] j ] (input data for 2R 11U ), we obtain the the input data. 
The output says 

2R llv x 2BT = 2xA2 + 2xA3-4xBl-4xB2 

+1 x CI + 1 x C2 + 2 x C3 + 0(h 2 ) , (7) 

where A2 = d^xhx^d^dyh^, A3 = d^dxhx^-d^dyh^, Bl = dudxh^-dxd^h^, B2 = 
d 2 hxu-dxd^h^ u ,Cl = d /Ji d v hxyd IJi d v h aa ,C2 = d 2 h flu -d 2 h^ lu ,C3 = d^d u h X x-d 2 h^ u , (Ta- 
ble 1, Appendix A). 

10 See O] for the present convention of gravitational quantities. 
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6 Discussion and Conclusion 



When we do the suffix contraction in usual algebraic softwares , such as REDUCE, 
we cannot help expressing the suffix permutation symmetry by explicitly writing up 
all possible permuted cases. For example, in REDUCE, Q defined in Sec. 2 is coded 
as 



INTEGER M,N; 
REAL Q; 

FOR ALL M,N LET DDH[M,N,M,N]=Q; 

FOR ALL M,N LET DDH[M,N,N,M]=Q; 



where an array DDH[M,N,L,S] expresses d^dphxa- It is very hard to generalize 
this approach to higher-power terms or to higher-rank tensors because it directly 
deals with all suffixes. The number of programming lines for the definition of each 
invariant increases rapidly. The present proposed approach do not treat suffixes 
directly, but the topology of graphs (, in other words, connectivity of suffix lines). 
The good efficiency of the algorithm based on the idea has been checked by some 
concrete programs. 

In this paper we have explained the new algorithm using a simple invariants 
of (ddh) 2 . Its generalization to other types can be done in a similar way. As the 
order of invariants increases, additional indices are required to discriminate between 
invariants (graphs). In we have used a generalized version of this program 
and have obtained the weak-field expansion of different products of three Riemann 
tensors. We have treated (<9<9/i) 3 -invariants (totally 90 terms) and the four indices 
are additionally required: 1) Vertex- Type Order, 2) Number of <9<9-vertices and h- 
vertices, 3) Number of crossing, 4) Disconnectivity. Most general SO(n) tensors 
(including invariants) constructed from dh, ddh, dddh, ■ ■ ■ can be treated similarly. 
The point in this algorithm is to find good indices to discriminate between invariants. 

The present algorithm is a new efficient way to do the suffix contraction using a 
computer. It can be applied not only to gravitational theories but also to general 
field theories as stated in Sec.l. One of interesting applications to physics is the Weyl 
anomaly calculation in gravitational theories. It is known that (ddh) m - type tensor 
calculation is sufficient for the Weyl anomaly calculation in 2m-dim gravitational 
theories. It was confirmed that Weyl anomaly in 4 dim quantum gravity is obtained 
by the present program [p~I . 
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♦*"*♦ * * * * i* " *jt 

Al A2 A3 

Figure 14: Graphs of (<9<9/i) 2 -invariants, loop no=l. 

Bl B2 B3 



: I . : I • 




B4 QQ 



Figure 15: Graphs of (<9<9/i) 2 -invariants, loop no=2. 

Appendix A. Graphical Representation of 
(<9<9/i) 2 -Invariants 

In the weak field expansion of n-dim Euclidean gravity : = S^ u + h^ u , \ h\ <C 1, 
we must generally treat global SO(n)-tensors which are composed of h^, d a h^ ui 
dadph^, ■ ■ -. Let us focus here on those tensors which are composed only of dadph^. 
We introduce its graphical representation as shown in Fig.l in Sec. 2 of the text. 
Contraction of suffixes is graphically represented by connecting the dotted lines 
with the same suffix. For example see Fig. 2 and 3 of Sec. 2 

In Sec. 2 we have graphically presented <9/i-invariants (Fig.3). They are P=d 2 h , 
Q—d^duh^ . Similarly we can list up all independent (<9<9/i) ^invariants as follows. 
They are grouped by the number of the suffix-loop: loop no =1, Fig. 14; loop no =2: 
Fig. 15; loop no =3: Fig. 16; loop no =4: Fig. 17;. 
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CI C2 



♦ ♦ ♦ * ♦ ♦ ♦ ♦ ♦ 



C3 PQ 

Figure 16: Graphs of (<9<9/i) 2 -invariants, loop no=3. 



• • - 



pp 

Figure 17: Graphs of (<9<9/i) 2 -invariants, loop no=4. 
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We have totally 13 invariants (10 connected, 3 disconnected). Their literal ex- 
pressions are listed in Table 1 of the text(Subsec.4.4). 
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Appendix B. Input and Output Data of Sample 
Calculations 

Example 1, Input Data 

4 
1 

12 3 4 
-1 

2 3 14 
-1 

14 3 2 
1 

4 3 12 

4 

1 

12 3 4 
-1 

2 3 14 
-1 

14 3 2 
1 

4 3 12 



Example 1, Output Data 
start of MAIN 

Input Check 
termnoa=4 
termO, weight =1 



d2ha[0 
d2ha[0 
d2ha[0 
d2ha[0 : 



terml, weight =-1 



d2ha[l 
d2ha[l 
d2ha[l 
d2ha[l 



d2ha[2 
d2ha[2 
d2ha[2 



term2, weight =-1 



0]=1 
1]=2 
0]=3 
1]=4 



0]=2 
1]=3 
0]=1 
1]=4 



0]=1 
1]=4 
0]=3 
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d2ha[2] [1][1]=2 
term3, weight =1 
d2ha[3][0][0]=4 
d2ha[3] [0][1]=3 
d2ha[3] [1][0]=1 
d2ha[3] [1][1]=2 
termnob=4 
termO, weight =1 
d2hb[0][0][0]=l 
d2hb[0] [0][1]=2 
d2hb[0] [1][0]=3 
d2hb[0] [1][1]=4 
terml, weight =-1 
d2hb[l][0][0]=2 
d2hb[l][0][l]=3 
d2hb[l] [1][0]=1 
d2hb[l][l][l]=4 
term2, weight =-1 
d2hb[2] [0][0]=1 
d2hb[2][0][l]=4 
d2hb[2][l][0]=3 
d2hb[2] [1][1]=2 
term3, weight =1 
d2hb[3] [0][0]=4 
d2hb[3] [0][1]=3 
d2hb[3][l][0]=l 
d2hb[3] [1][1]=2 



FINAL ANSWER 

-8 Al 

A2 

A3 

Bl 

B2 

4 B3 

4 B4 

QQ 

CI 

C2 

C3 

PQ 

PP 



Example 2, Input Data 



4 
1 

12 3 3 
-1 

13 3 2 
-1 

2 3 3 1 
1 

3 3 12 
4 

1 

12 4 4 
-1 

14 4 2 
-1 

2 4 4 1 
1 

4 4 12 



Example 2, Output Data 
start of MAIN 



Input Check 




termnoa= 


4 




termO,wei 


ght: 


=1 


d2ha[0] [0] 


[0] = 


--1 


d2ha[0] [0] 


[1] = 


-2 


d2ha[0] [1] 


[0] = 


-3 


d2ha[0] [1] 


[1]= 


-3 


terml,wei 


ght: 


=-1 


d2ha[l][0] 


[0] = 


--1 


d2ha[l] [0] 


[1] = 


-3 


d2ha[l] [1] 


[0] = 


-3 


d2ha[l][l] 


[1]= 


--2 


term2,wei 


ght: 


=-1 


d2ha[2] [0] 


[0] = 


-2 


d2ha[2] [0] 


[1]= 


-3 


d2ha[2] [1] 


[0] = 


-3 


d2ha[2] [1] 


[1]= 


-1 


term3, weight - 


= 1 



Aoh o tqi mi 
aziia[oj [UJ 


mi — 

UJ- 


-o 


H9WVI D 


Ml- 


-o 


U-Zrlldl OJ 111 


rni- 

UJ- 


-1 

- 1 


H9hfl fulfil 

U-ZdldlOJ 1 ±1 


[11- 


-9 


tprrnnnn — 

Lid \ J — 


4 




tt;IIIlU,Wt;l 


ght 


_1 
— 1 


ri9hhfni rnl 

U.Z,11U 1 UJ 1 UJ 


mi — 

[UJ- 


-1 

- 1 


H2hhfni rnl 


[11- 


-2 




rni — 


A 


d2hb 1 


[11- 


A 


terml,wei 


ght 


— i 


HPhhTlirnl 


mi — 

[UJ- 


-i 

- J. 


H2hhfl rn 


[11- 


-4 


LlZdlU 111 111 


rm- 

[UJ- 


-4 


d2hb 1 1 


[il- 


-0 

- Zj 


term2,wei 


ght 


— 1 


d2hb[2] [0] 


[0] = 


-9 


d2hb[2] [0] 


[1] = 


-4 


d2hb[2][l] 


[0] = 


-4 

-4: 


d2hb[2][l] 


[1]= 


-i 

- J. 


term3, weight 


=1 


d2hb[3] [0] 


[0] = 


A 


d2hb[3] [0] 


[1]= 


A 


d2hb[3] [1] 


[0] = 


A 


d2hb[3] [1] 


[1]= 


=2 



FINAL ANSWER 

Al 

2 A2 

2 A3 

A Bl 

-4 B2 

B3 

B4 

QQ 

1 Cl 

1 C2 

2 C3 
PQ 
PP 
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Figure Captions 

• Fig.l 4-tensor d^d v h a p 

• Fig.2 2-tensors of (a) d 2 h a/ 3 , (b) d^d u h aa and (c) d^dph a p 

• Fig.3 Invariants of P = d^d^h aa and Q = d a dph a p . 

• Fig. 4 4-tensor of did 2 h 3 4. 

• Fig. 5 Invariants of P = d\dih 22 and Q = d\d 2 h\ 2 ■ 

• Fig. 6 8-tensor of did 2 h 34 ■ d 5 d 6 h 7S . 

• Fig. 7 4-tensor of did 2 h 34: • d 5 d 6 h 3 4. 

• Fig. 8 Invariant of 53 = d\d 2 h^ ■ did 2 h 34: . 

• Fig.9 loopstream j ][ ][ ] for B3. 

• Fig. 10 Invariant of Bl. # of vertices of 0-th and 1-th loop are 1 and 3 respec- 
tively. 0-th loop is a tadpole. 

• Fig. 11 tadpoleno and tad[ ][ ] for CI, C2 and C3. 

• Fig. 12 Change of i (bond number) and j (vertex-type number). 
Arrows indicate directions of tracings. 

• Fig. 13 absdeli[ ] and absdelj [ ] for Al, A2 and A3. 

• Fig. 14 Graphs of (<9<9/i) ^invariants, loop no=l. 

• Fig. 15 Graphs of (<9<9/i) ^invariants, loop no=2. 

• Fig. 16 Graphs of (<9<9/i) ^invariants, loop no=3. 

• Fig. 17 Graphs of (<9<9/i) ^invariants, loop no=4. 
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